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In a boost-invariant hydrodynamic model , we have analyzed <j> meson production in Pb+Pb and 
Au+Au collisions, in the centre of mass energy range yfs= 6-200 GeV. Hydrodynamic evolution 
is governed by a lattice based equation of state with a confinement-deconfinement cross over at 
T co =196 MeV. We also look for the threshold energy above which the confined matter undergoes a 
deconfinement transition. Analysis indicate that above a threshold energy, y / s th = 13.65± 3.06 GeV, 
the fluid medium crosses over from a confined phase to a deconfined phase. 
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I. INTRODUCTION 

Lattice simulations 0, 0] of Quantum Chromody- 
namics predict that nuclear matter, at extreme condi- 
tions, e.g. high temperature or pressure, can undergo 
a confinement-deconfinement cross-over transition. The 
deconfined matter, comprising quarks and gluons, even 
at very high temperature, is not totally free from interac- 
tion and one calls the matter strongly interacting Quark 
Gluon Plasma (sQGP). Heavy ion collisions at relativis- 
tic energies provides a convenient method to produce 
and study properties QCD matter at extreme conditions. 
Recent experiments in Au+Au collisions at RHIC 
produced convincing evidences that in central and mid- 
central Au+Au collisions, a hot, dense, strongly inter- 
acting, collective QCD matter is created. Whether the 
matter can be characterized as the lattice QCD [l|, Q 
predicted Quark-Gluon-Plasma (QGP) or not, is still a 
question of debate. 

For long, strangeness enhancement is considered as a 
signature of QGP formation [7]. In QGP environment, 
gg — > ss is abundant. If not annihilated before hadro- 
nisation, early produced strange and anti-strange quarks 
will coalesce in to strange hadrons and compared to ele- 
mentary pp collisions, strange particle production will be 
enhanced. However, strangeness enhancement could also 
be obtained in apurely hadronic scenario, mainly due to 
'volume effect' [8Hlll|. Strangeness production in small 
volume elementary pp collisions can be 'canonically' sup- 
pressed due to 'strict' strangeness conservation [a-lllj. In 
bigger volume AA collisions, locally, strangeness conser- 
vation condition can be relaxed to produce strange par- 
ticles. In the language of statistical mechanics, while 
canonical ensemble is applicable in pp collision, grand 
canonical ensemble is more appropriate in heavy ion col- 
lisions. Additionally, strange particle phase space ap- 
pears to be undersaturated in elementary pp or periph- 
eral heavy ion collisions [H, [Hj]. The deviation from 
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strangeness equilibrium is parameterized by strangeness 
undersaturation factor 73. In a statistical model, fits to 
particle multiplicities are unsatisfactory, unless 7s is not 
accounted for. 75 shows an increasing trend from AGS 
to RHIC energy. At RHIC, 75 ~ 1. A possible interpre- 
tation is that the total amount of strangeness available 
at the hadronisation is determined in the pre-hadronic 
stage. The change in 75 between pp and A+A collisions 
then reflect the difference between initial conditions of 
the respective systems. 

Recently STAR collaboration published their measure- 
ments of <f> production in Au+Au and Cu+Cu collisions 
at RHIC energy (yi=200 GeV) dM^. Both in Au+Au 
and Cu+Cu collisions, (f> meson production is enhanced. 
The enhancement factor is ~ 2. 4>(ss) mesons are hid- 
den strange particles, unaffected by the 'canonical' sup- 
pression. STAR data on (f> meson production is consis- 
tent with models with recombination of thermal strange 
quarks, indicating that in Au+Au collisions a robust par- 
tonic system in created [17l |. Interestingly, strangeness 
enhancement is observed even at much lower energy. At 
SPS, NA49 collaboration measured 4> meson production 
in 20A, 30A, 40A, 80A and 158A GeV Pb+Pb collisions 
|18| . 4> meson production is enhanced also at SPS en- 
ergy. Enhancement factor, ~3-4, is even larger than 
that in RHIC collisions. One naturally wonders about 
the sources of <fi meson enhancement in low SPS en- 
ergy collisions. Is it due to formation of a deconfined 
medium like in Au+Au collisions at RHIC? NA49 col- 
laboration analyzed the data in the statistical hadronisa- 
tion model 12]. Statistical hadronisation model, includ- 
ing the strangeness undersaturation factor 7s, provides a 
good description of the data. Apparently, as in Au+Au 
collisions at RHIC, initial partonic content of the system 
also drive the <p meson production in Pb+Pb collisions 
at SPS. 

Statistical hadronisation model is not a dynamic 
model. A statistical model can not comment whether or 
not the initial system is in a deconfined state. Relativistic 
hydrodynamics is the only dynamical model, which can, 
within certain limitations, comment on the initial condi- 
tion of the fluid. It is assumed that in the collision a fire- 
ball is produced. Constituents of the fireball collide fre- 
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quently to establish local thermal equilibrium sufficiently 
fast and after a certain (proper) time T;, hydrodynam- 
ics become applicable. If the macroscopic properties of 
the fluid e.g. energy density pressure, velocity etc. are 
known at the equilibration time Tj , the relativistic hydro- 
dynamic equations (i.e. energy-momentum conservation 
equations) can be solved to give the space-time evolu- 
tion of the fireball till a given freeze-out condition such 
that interactions between the constituents are too weak 
to continue the evolution. Using suitable algorithm (e.g. 
Cooper-Frye) information at the freeze-out can be con- 
verted into particle spectra and can be directly compared 
with the experimental data. Thus, hydrodynamics, in an 
indirect way, can characterize the initial condition of the 
medium produced in heavy ion collisions. Hydrodynam- 
ics equations are closed only with an equation of state 
(EOS) and one can investigate the possibility of phase 
transition in the medium. 

In the present paper, in an ideal hydrodynamic model, 
we have analyzed the NA49 data on 4> meson produc- 
tion in 20-158 A GeV central Pb+Pb collisions and the 
STAR data [ll-OJ] on cj> meson production in y/s=62, 130 
and 200 GeV central Au+Au collisions, effectively scan- 
ning the cm. energy range -^5=6-200 GeV. An impor- 
tant ingredient of the model is a lattice based EOS with 
a confinement-deconfinement cross over transition at 
T c0 =196 MeV. Our analysis indicate that ideal hydrody- 
namics, with the lattice based EOS, reasonably well ex- 
plain the 4> meson production in the energy range, a/s=6- 
200 GeV. Analysis also indicate that above a threshold 
energy yfs th = 13.65 ± 3.06 GeV, the medium crosses 
over from a confined phase to a deconfined phase. It was 
also indicated that the transition is rather sharp, occur- 
ing over a narrow energy range Ay/s th = 11.28 ± 5.93 
GeV. 

The paper is organized as follows: in section [TTl we 
briefly describe the construction of a lattice based equa- 
tion of state. Initial conditions used to compute fluid evo- 
lution etc., are also discussed in section [TTl Experimental 
data e.g. <j> meson's multiplicity, py-spectra, mean px 
are analyzed in section Mil Summary and conclusions 
are given in section IIVI 
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FIG. 1: Black circles are lattice simulation 0] for entropy 
density. The dash line is a parameterisation of lattice simu- 
lation by Eq[T] The solid line is the entropy density in the 
lattice+HRG model (see text). 



phase transition. However, lattice simulations |2| indi- 
cate that confinement to deconfinement transition is a 
cross over, rather than a 1st or 2nd order phase transi- 
tion. It is then essential that hydrodynamic simulations 
are done with EOS with cross-over transition rather than 
with EOS with 1st or 2nd order transition. Huovinen 

21] reported an 'ideal' hydrodynamic simulation with a 
cross-over phase transition. He concluded that an EOS 
with 1st order phase transition better explain the exper- 
imental data (e.g. elliptic flow of proton or antiproton) 
than an EOS with cross over transition. 

Presently, we have used an EOS based on recent lattice 
simulation of Cheng et al. 2]. In Fig[IJ simulation results 

2] for the entropy density is shown. The dotted line in 
FigQ]is a parameterisation of the entropy density. 



± E = a + [P + 1 T][l + tanh I ^±], (1) 



II. HYDRO DYNAMICAL EQUATIONS, 
EQUATION OF STATE AND INITIAL 
CONDITIONS 



with a=0.64, 0=6.93, 7=0.55, T c =196 MeV, AT = 
0.1T C . From the parametric form of the entropy den- 
sity, pressure and energy density can be obtained using 
the thermodynamic relations, 



Assuming boost-invariance, we have solved the energy- 
momentum conservation equation, dfjT^ v = in (r = 
y/t 2 — z 2 , x, y, r/ s — \ In jzr^) coordinates, using the code 
'"AZHYDRO", details of which can be found in [H. As 
mentioned earlier, we have used a lattice based equation 
of state. Equation of state (EOS) is one of the most im- 
portant inputs of a hydrodynamic model. Through this 
input macroscopic hydrodynamic models make contact 
with the microscopic world. Most of the hydrodynami- 
cal calculations are performed with EOS with a 1st order 



p(T) = [ s(T)dT (2) 
Jo 

e{T) = Ts-p. (3) 

We complement the lattice simulated EOS by a 
hadronic resonance gas (HRG) EOS comprising all the 
resonances below mass 2.5 GeV. In Fig(TJ the solid line 
is the entropy density of the '"lattice +HRG" ' EOS. The 
entropy density of the complete EOS is obtained as, 
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FIG. 3: Ratio of the mean <j> meson multiplicity per par- 
ticipant over cube of the initial temperature, normalised by 
the participant number (N part ) and rapidity width (AY), 
in Pb+Pb and Au+Au collisions as a function of collisions 
energy. The solid line is a fit to the ratio with the form, 



all + tanh 



FIG. 2: In panels (a)-(e), transverse momentum spectra for 
cj> mesons in Pb+Pb collisions at 7.6, 8.8, 12.3 and 

17.3 GeV are shown. STAR measurements of pr spectra for 
4> in y / s=62,130 and 200 GeV Au+Au collisions are shown 
in panels (f)-(h). The lines are respectively are ideal hydro- 
dynamic fit to the <j> spectra. The blue line in panel (a), is 
a fit obtained to the data when confinement-deconfinement 
transition occur at T co = 160 MeV. 



s(h,x,y) = e [0.75N part (b,x,y) + 0.25iV coH (b,ir,y)], 



where b is the impact parameter of the collision. N, 



(5) 



part 



and N co u are the transverse profile for the average partic 
ipant and collision number respectively. N part and N co u 
are calculated in a Glauber model, with nuclear density. 



s = 0.5[1 + tanh(x)]sHRG + 0.5[1 — tanh(x)]s lattice 

(4) 

with x = T A ^ C . Compared to lattice simulation, entropy 
density in HRG drops slowly at low temperature. It is 
consistent with the observation that at low temperature, 
trace anomaly, g ~4 P drops faster in lattice simulation 
than in a HRG model [2] . It is difficult to resolve whether 
the discrepancy between lattice simulations at low tem- 
perature and HRG model is due to failure of HRG model 
at lower temperature or due to the difficulty in resolving 
low energy hadron spectrum on a rather coarse lattice 

We note that, lattice simulations [2[ were performed 
with zero baryon density /z_b=0. /is=0 is only approx- 
imately valid in RHIC energy collisions. In lower SPS 
energy collisions, there can be considerable stopping and 
Hb can be appreciably different from zero. However, lat- 
tice simulations with finite baryon density have yet not 
reached the state of application. We therefore continue to 
use the Hb=® lattice based EOS, even at low SPS energy 
collisions. Later, we will comment on the possible effects 
finite baryon density EOS can have on our analysis. 

We assume that at the initial time Ti=0.6 fm, the initial 
energy density of the fluid is distributed as [20[ 



P^) = L + e (r-fl)/„ ' ( 6 ) 

with R=6.624 (6.37) fm and a=0.549 (0.54) fm for Pb 
(Au) nucleus. p$ is determined from the normaliza- 
tion condition, / p(r)d 3 r = A (208 for Pb, 197 for Au). 
Glauber model calculation also require the total inelastic 
NN cross section <Ji n ei- NN cross-section varies slowly 
with energy. Throughout the energy range, t/s=6.3-17 
GeV, we assume cr ine i— 32 mb. For -^5=62, 130 and 200 
GeV Au+Au collisions, we have used, <x;, ie z=36, 40 and 
42 mb respectively. We also assume that at the initial 
time Tj, the transverse fluid velocity is zero. £n m Eqf5] 
corresponds to central energy density in zero impact pa- 
rameter collisions. We fix it to reproduce the experimen- 
tal pt spectra of <f> mesons from a freeze-out surface at 
T F =150 MeV [H. 

III. RESULTS 

Details of the data analysed are given in tableHl NA49 
collaboration measured transverse momentum spectra 
(dP^) of <t> mesons in ^5=6.3, 7.6, 8.8, 12.3 and 17.3 
GeV Pb+Pb collisions. The measured spectra are shown 



4 



TABLE I: Details of the data analysed are given. The Pb+Pb data are from the NA49 experiment [l8| . The Au+Au data are 
from the STAR experiment [l5|]. b is the impact parameter corresponding to the centrality cut. The fitted values of the central 
energy density, eo and corresponding /N of the fit, ratio of experimental cj> multiplicity over the hydrodynamic predictions 
< TV" > ex I < N > Tir and the ratio of experimental meson mean pr over the predicted mean pr, < pt > ex / < Pt > th are 
shown. 



System 


(GeV,) 


collision 
Centrality 


rapidity 
acceptance 


b 

ffrr,\ 

(imj 


eo 

(LreV/jm ) 


X 2 /N 


Q( X 2 /N,N) 


< iv > 
<N>TH 


<PT> TH 


Pb+Pb 


6.3 


0-7.2 % 


0-1.8 


2.71 


1.32 


0.69 


0.98 


2.39 ± 1.43 


1.22 + 0.09 


Pb+Pb 


7.6 


0-7.2 % 


0-1.8 


2.71 


1.22 


3.97 


0.55 


2.00 + 0.60 


1.73 + 0.12 


Pb+Pb 


8.8 


0-7.2 % 


0-1.5 


2.71 


1.77 


4.52 


0.92 


1.83 + 0.33 


1.33 + 0.03 


Pb+Pb 


12.3 


0-7.2 % 


0-1.7 


2.71 


2.87 


5.08 


0.75 


1.61 ± 0.26 


1.32 + 0.04 


Pb+Pb 


17.3 


0-5.0 % 


0-1.0 


2.3 


7.15 


9.23 


0.51 


1.09 + 0.05 


1.29 + 0.03 


Au+Au 


62.0 


0-20.0% 


-0.5-0.5 


4.43 


18.16 


1.52 


0.99 


1.17 + 0.22 


0.93 + 0.06 


Au+Au 


130.0 


0-11.0% 


-0.5-0.5 


3.2 


21.62 


13.94 


0.12 


1.50 + 0.21 


1.06 + 0.06 


Au+Au 


200.0 


0-5.0 % 


-0.5-0.5 


2.3 


31.2 


1.02 


0.99 


1.03 + 0.09 


0.91 + 0.06 



in five panels (a)-(e) of FigH In Vs=6.3-12.3 GeV 
Pb+Pb collisions, spectra are measured in 0-7.2% cen- 
trality collisions, while in ^=17.3 GeV Pb+Pb colli- 
sions centrality cut was 0-5%. STAR measurements for 
(j> mesons spectra in y/s—62, 130 and 200 GeV Au+Au 
collisions are shown in panels (f), (g) and (h). Central- 
ity cuts for y/s=62, 130 and 200 GeV Au+Au collisions 
are 0-20%, 0-11% and 0-5% respectively and the rapidity 
acceptance is \y\ < 0.5 . 

Before we proceed further, it is important to mention 
that in SPS energy Pb+Pb collisions, the assumption of 
boost-invariance is not valid. NA49 collaboration [l8[ 
measured rapidity distribution of <j) mesons. In the en- 
ergy range, -^=6.3-17.3 GeV, the rapidity distribution 
do not show plateau like structure, rather the distribution 
could be fitted with a Gaussian, width of the Gaussian 
increasing with energy 18]. 3+1 dimensional hydrody- 
namical model is then required to analyze NA49 data. 
The boost-invariant hydrodynamics is used in the spirit 
of local density approximation. The underlying assump- 
tion is that around each rapidity y, a narrow window 
exist where dN/dy is approximately constant. As will 
be shown later, with progressive energy, the assumption 
become more and more accurate. 

The NA49 and STAR measurements for the <fi meson 
Pt spectra arc fitted in the model by varying the central 
energy density e - In Fig HI the solid lines are the best fit 
obtained to the data. In table|H best fitted values of the 
central energy density and corresponding x 2 /N , for the 
systems analysed, are noted. In the x 2 calculations, we 
have included the statistical and systematic errors [22l |. 
To test the goodness of fit, in tableJH we have noted the 
probability Q( X 2 /N,N), 

PX — OO 

Q((x 2 /N) mm ,N)= X 2 (x,N)dx, (7) 

Jx=( X 2 /N) min 

that a random sampling from a Normal distribution 
would exceed (x 2 /N) m i n . A statistically significant fit 



require Q > 0.5. With the exception of <fi meson data 
in "^5=130 GeV Au+Au collisions (quality of which is 
poor), all the data sets satisfy the requirement. Hydro- 
dynamic model, with the lattice+HRG EOS, reasonably 
well explain <fi meson spectra over a wide range of energy, 
v^=6-200 GeV. 

In tablelH the ratio of the experimental <p meson's 
peak multiplicity over the predicted multiplicity ( ) 
, as well as the ratio of the experimental mean pt over 
the predicted mean px ( fepT? ) > are shown. The er- 
rors in (n}th and ^^th are experimental only. Dis- 
agreement with model predictions for (f) multiplicity and 
experiments increases as the collision energy is lowered 
(with the exception of 130 GeV Au+Au data). Con- 
tinual increase of discrepancy between the model predic- 
tions and experiments with lowering energy indicate that 
the assumption of boost-invariance is violated more and 
more at lower energy. In the energy range y / s=6-12 GeV 
Pb+Pb collisions, NA49 measured 4> mesons in the ra- 
pidity range Ay = — 1.5(1.8). We have fitted the px 
spectra integrated over the rapidity range. Over this ra- 
pidity range, <fi meson multiplicity changes considerably. 
Naturally, evolution of the fluid, which fit (rapidity) in- 
tegrated pt spectra, underpredict multiplicity at central 
rapidity. <fi mesons mean px at SPS energy is also under- 
predicted in the model. 

If beyond a threshold energy, there is a cross over tran- 
sition from the confined hadronic matter to the decon- 
fined QGP, there will be a rapid increase in the degener- 
acy of the medium (e.g. see Fig|T]). Will there be a signa- 
ture of the transition in the energy range -^=6-200 GeV 
scanned by the 4> mesons? Ratio of the entropy density 
over cube of the temperature is a measure of the degener- 
acy of the medium. </> meson's multiplicity is expected to 
be proportional to the initial entropy density. The ratio 
of <j) multiplicity over the cube of the initial tempera- 
ture of the fluid is then proportional to the degeneracy 
of the medium. In Fig. 3, collision energy dependence of 
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the ratio of the experimental <fi meson multiplicity over 
the cube of the spatially 'averaged' initial temperature 
< N > /[.5 ■ N part ■ AY ■ T 3 } is shown. Multiplicities are 
normalised by the participant numbers (N part ) and ra- 
pidity gap (AF), to account for the system size, different 
centrality bins and different phase space in Au+Au and 
Pb+Pb collisions. The ratio is approximately constant 
at RHIC energy v / s=62-200 GeV. The ratio is approxi- 
mately constant also at the low SPS energy, y/s < 12.3 
GeV. From low SPS energy to RHIC energy, the ratio 
increases by a factor of ~ 4. It may be noted (see Figfl]) 
that in the lattice+HRG EOS, entropy over cube of the 
temperature also increases by a similar factor (~ 4.5) in 
the temperature range 150-400 MeV. 

In FiglJl the solid line is a fit obtained to the ra- 
tio with an analytical form for the step function, a[l + 

tanh V T~/' th }. Fitted values are, a = 1.05 ± 0.21, 

^^=13.65 + 3.06 GeV and Ay/s th = 11. 28 ±5.93 GeV. 
The threshold energy \/s th can be determined within ~ 
20% accuracy. Presently no experimental data exists in 
the energy range 17.3-62 GeV, between the top of the 
SPS energy and bottom of the RHIC energy. Uncertainty 
in the threshold energy could be reduced if the gap is 
filled. Future STAR energy scan program at RHIC may 
help to narrow down the threshold energy [23|, [24| ■ It is 
also heartening to find the threshold energy close to the 
top of the SPS energy, where first evidence of QGP for- 
mation was obtained [251 ] . The transition is also rather 
sharp, over a narrow energy range Ay/s th ps 12 GeV, the 
confinement-deconfincment transition occur. If the gap 
between \/s=l7.3-62 GeV is filled in future RHIC energy 
scan programme, width of the transition may be further 
reduced. 

In the present analysis we have assumed a baryon free 
fluid and used the lattice based (J-b—0 equation of state. 
While at RHIC energy collisions, is approximately 

valid, at SPS energy collisions, /is ^ 0. In a finite baryon 
density fluid, the confinement-deconfincment transition 
temperature decreases [26j-|28[ . Also there is a possibility 
that at some hb, the confinement-deconfinement transi- 
tion become a 1st order phase transition, with a critical 
point at some intermediate the baryon density. At the 
critical point, conserved charges (e.g. baryon density) 



will have large fluctuations 29]. The effect of finite hb 
on <f) meson production will be an interesting study, which 
we could not do presently. However, following comments 
can be made. Other conditions remaining unchanged, in 
a finite [Lb fluid, <fr yield will decrease as the cross-over 
temperature decreases. Then to reproduce the experi- 
mental multiplicity, initial temperature of \ib fluid 
will be higher than fis = fluid. As an example, we have 
fitted 4> meson px spectra in -^=6.3 GeV Pb+Pb col- 
lisions with cross-over transition at T co =160 MeV. The 
best fit to the data is obtained with central energy den- 
sity Eq—IA GeV/ fm 3 , ~ 6% higher than the value re- 
quired for X" co =196 MeV. The fit is shown as the blue 
line in Figj2k. It can not be distinguished from the fit 
obtained when T co =196 MeV. With reduced initial tem- 
perature in SPS energy collisions, the ratio < N > /T 3 
will decreases, and the confinement-deconfinement tran- 
sition will be sharper than estimated presently. 



IV. SUMMARY AND CONCLUSIONS 

To summarise, assuming that in high energy heavy 
ion collisions, a fluid medium is created, which become 
amenable to hydrodynamic description after Ti=0.6 fm, 
we have analysed the recent NA49 data on <j> meson pro- 
duction in Pb+Pb collisions in the energy range \/s=6.3- 
17.3 GeV and the STAR data on <\> meson production 
in ^=62,130 and 200 GeV Au+Au collisions. The 
hydrodynamic evolution is governed by a lattice+HRG 
EOS, with a confinement-deconfinement cross-over tran- 
sition at T co =196 MeV. While the confined phase of the 
EOS is modelled by a non-interacting hadronic resonance 
gas, the deconfined part is modelled after a recent lat- 
tice simulation. (f> meson production in nuclear colli- 
sions, over the energy range i/s=6. 2-200 GeV, are rea- 
sonably well explained in the model. We have also tried 
to extract the threshold energy above which a confined 
medium crosses over to a deconfined medium. Energy 
scan by the <f> meson indicate that above a threshold en- 
ergy \fs th = 11.28 ± 5.93 GeV, the medium crosses over 
from a confined to a deconfined phase. 
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